The Effect of Symmetry Lowering on the Dielectric Response of BaZrOs 
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We use first-principles density functional theory calculations to investigate the dielectric response 
of BaZrOs perovskite. A previous study [Arkbarzadeh et al. Phys. Rev. B 72, 205104 (2005)] 
reported a disagreement between experimental and theoretical low temperature dielectric constant 
e for the high symmetry BaZrOa structure. We show that a fully relaxed 40-atom BaZrOs structure 
exhibits Oe octahedral tilting, and e that agrees with experiment. The change in e from high- 
symmetry to low-symmetry structure is due to increased phonon frequencies as well as decreased 
mode effective charges. 



I. INTRODUCTION 

Dielectric materials are important for wireless commu- 
nications technology. These devices require a high di- 
electric constant, e, and low dielectric loss fl|. Barium 
zirconate, BaZrOa (BZ), is one of the constituent ma- 
terials in the electroceramic capacitors used in wireless 
communications. BZ is a classic ABO3 perovskite dielec- 
tric material that is both chemically and mechanically 
stable. 

According to the well-established tolerance factor ar- 
gument, BZ should have a stable cubic structure. Toler- 
ance factor t is given by 
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where Ra-o is the sum of A and O ionic radii and Rb-o 
is the sum of B and O ionic radii. Tolerance factor t < 1 
usually leads to the rotation and expansion of the B- 
Og octahedra. Such octahedral rotations often generate 
a low temperature anti-ferroelectric (AFE) phase (e. g. 
PbZrOs). If i > 1, the B-Oe octahedra are stretched from 
their preferred B-0 bond lengths, promoting B-cation 
distortions by creating room for the B-cations to move 
off-center. Therefore, simple perovskites with i > 1 are 
usually ferroelectric (FE). For BZ, the sizes of the Ba and 
Zr ions (i?Ba=l-61 A and i?zr=0.72 A), in their O12 and 
Og cage (i?o = l-35 A) exactly balance, leading to t = 1, 
so the cubic structure should have no driving force to 
deform to a lower symmetry. It also has a stable dielectric 
constant over a wide range of temperatures 0| • 

The total dielectric constant is a sum of two contribu- 
tions, related to the electronic (eoo) and ionic polarizabil- 
ities (e^i). For most perovskites, Coo is small (« 5) and 

is the dominant contribution at low frequency. Recent 
experimental studies of BZ at low temperature show that 
e is 47 at K and some interesting quantum effects are 
present, though not large in magnitude Q. Reference 3 
also presents a DFT study which calculates e of BZ to 
be 65 at K, overestimating experimental results. DFT 
inaccuracies have been known for some time in describ- 
ing ferroelectric systems Q . It is however surprising that 



DFT is inaccurate for a system as simple as paraelectric 
BZ. 

Here, we show that at low temperatures, the ground 
state structure of BZ is not accurately represented by a 
five atom cell. A larger supercell with lower symmetry is 
energetically favored. The calculated value of e for this 
low symmetry structure is in agreement with experimen- 
tal results. Most importantly, the relationship between 
the breaking of structural symmetry and the change in 
dielectric response is elucidated. 



II. 



METHODOLOGY 



In this study, two first principles codes are used. An 
in-house solid state DFT code used in previous stud- 
ies H, Q and the ABINIT software package are used 
to relax the ionic positions and lattice constants. Local 
density approximation (LDA) of the exchange correlation 
functional and a 4 x 4 x 4 Monkhorst-Pack sampling of 
the Brillouin zone Q are used for all calculations. All 
atoms are represented by norm-conserving optimized 
designed nonlocal pseudopotentials. All pseudopo- 
tentials generated using the OPIUM code ll]. The cal- 
culations are performed with a plane wave cutoff of 50 Ry. 

Once the structure is fully relaxed, response func- 
tion 12, 13] calculations are performed with ABINIT to 
generate Dap{i,j), the mass weighted dynamical matrix. 
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Once normalized, each eigenvalue is v'^, the frequency 
squared of a normalized eigenvector a^. Born effective 
charge tensors, Z*^p are also calculated for each atom. 
Each mode /x has a mode effective charge, Z*^, defined 
as 
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Contributions to the dielectric response arise only from 
the IR active modes (l^. These modes have a non-zero 
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TABLE I: Forty-atom BZ phonon data for IR active modes in 
the high-symmetry cubic structure. Frequencies are in cm~^, 
mode eflective charges, Z*, are in electrons. The contribu- 
tions to the e are calculated in SI units. Also reported are the 
motions which create the IR active mode. 
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and are used to calculate e^, the contribution to the 
dielectric constant from mode n |15| as 



(4) 



where the total ionic contribution is 



(5) 



III. 



RESULTS AND DISCUSSION 



A. Five-Atom BZ calculations 

The ground state five atom BaZrOa is a cubic per- 
ovskite structure where all ions occupy high symmetry 
positions. The relaxed DFT lattice constant of this cell 
was a =4.157 A. This structure was used to calculate 
both a phonon band structure (Figure and e. The 
high-symmetry BaZrOa structure contains three IR ac- 
tive phonon modes in each direction. 

The first set of phonon frequencies, occurring at 96 
cm"^ contribute the most to e, as shown in Table ^ 
These three modes are mainly large Ba displacements 
against its O12 cage into which some Zr displacement 
opposite its Og cage is mixed (Last mode). The ratio of 
Ba displacement to Zr displacement is about 9:1. How- 
ever, due to the much larger Z* of Zr (6.3 for Zr versus 
2.7 for Ba) even such a small Zr off-center displacement 
makes a significant contribution to the effective charge of 
the mode. The IR active modes are sometimes discussed 
in terms of purely A-0 or B-0 displacements, but that 
simple cation motion is not justified when describing the 
lowest frequency phonon modes of BZ. 

The second set of phonon frequencies occurs at 193 
cm~^, consisting of Zr motion (Slater mode). A small 
amount of Ba motion is also present, but the ratio of Zr 
off-centering to Ba off-centering in this mode is about 
20:1. Combined with the smaller Z* of Ba, this means 
that this mode is well-characterized as pure B-0. The 
final set of phonons, at 513 cm~^ are caused by asym- 
metric Og cage motions. These contribute the least to 
e. 



FIG. 1: The calculated phonon band structure of the five- 
atom high-symmetry BZ structure. Our results match well 
with the previously published BZ phonon band structure in 
reference 0]. In particular, our unstable mode at R is -79 
cm~^, in close agreement with their -63 cm~^ 



To obtain e, the contributions of from Table ^ are 
added to e^o, which was calculated as 4.94. The total 
zero temperature dielectric constant for the high symme- 
trystructure is 65. This value is almost double previous 
experimental investigations 2], but agrees well with re- 
cent theoretical investigations 3]. 

The disagreement between experimental and theoreti- 
cal e results suggests either an inaccuracy of LDA or an 
inaccurate representation of ground-state BaZrOa struc- 
ture by a five-atom high-symmetry unit cell. In Refer- 
ence Q , the authors present a phonon band structure of 
BZ, which displayed a soft mode at R, but none at F. 
Our phonon band structure of five atom BZ, presented 
in Figure n confirms that there is no ferroelectric soften- 
ing at F, but a rather large softening at R. The presence 
of an unstable mode indicates a driving force for symme- 
try breaking. To explore the effect of this, a larger BZ 
supercell is needed for accurate calculations. 



B. Forty-Atom BZ calculations 

The soft mode at R observed in the phonon band struc- 
ture was induced in the high symmetry structure, and the 
supercell was then relaxed. The relaxed low symmetry 
structure had lattice constants a = 8.306A, h = 8.313A, 
and c — 8.290A. Thus, the relaxed volume is only 0.39% 
less than the high symmetry structure. None of the 
cations shift from their high symmetry positions, consis- 
tent with an Og tilt mode. Octahedral tilt angles range 
between 3.6 and 4.2°. A subsequent phonon response 
calculation on the tilted pseudo-cubic structure showed 
that no negative frequency modes were present. Unlike 
our calculations, low temperature x-ray diffraction and 
neutron scattering experiments are unable to detect this 
tilting in BZ Q. Current experimental evidence is un- 
clear. The relevant X-ray reflections for oxygen-based 
symmetry breaking would be weak, due to the scatter- 
ing factor of O. We propose the tilted structure based 
on quantum mechanical energy minimization and conse- 
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TABLE II: Forty-atom low-symmetry BZ phonon data for 
IR active modes. Frequencies are in cm~^, mode effective 
charges, Z* , are in e. The contributions to e are calculated 
in SI units. For each IR active mode, the atomic motions are 
described, mode. 
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FIG. 2: Gamma-point phonon DOS of the high-symmetry 
BaZrOa structure, solid black line, and the fully relaxed 
pseudo-cubic supercell, dashed blue line. The most important 
changes are the upshift in Ba-O mode frequencies (from 96 
cm~^ to 100-130 cm~^) and the appearance of new modes 
around 300 cm^^. 



quent dielectric response (see below) and we suggest fur- 
ther experimental investigation to directly confirm the 
predicted structure. 

Analysis of the low-symmetry (Table ^ and Fig- 
ure |2Jl phonon modes show that there are more IR-active 
phonon modes than in the high-symmetry cell. The IR- 
active modes are no longer of pure x, y or z character. 
However, they can still be distinguished as Ba-O (Last 
mode), Zr-0 (Slater mode) or Og based. Likewise, the 
frequency ranking of the modes is still lowest for Ba-O, 
higher for Zr-0 and highest for the Og. The Ba-O modes 
are clustered around 115 cm~^, about 10% higher in fre- 
quency than for the high-symmetry cell. The energies 
of the Zr-0 and the Og modes are essentially unchanged 
from their respective 5-atom cell values. Interestingly, a 
new band of modes is generated by this structure around 
300 cm~^. do not contribute as much to as the three 
previously mentioned sets of modes. 

Comparison of the values in Tables I and II shows 
that the contribution of the mainly Ba modes to e has de- 
creased by half, whereas the contribution of the Zr modes 
and the Og modes has remained the same. The new 
modes around 300 cm~^ do not contribute significantly. 
The total is 45, in contrast with 60 calculated for the 
5-atom cell. We find that is still 4.94, and the total 
e is 50, in agreement with experimental data. Thus, we 
find that it is necessary to include the Og tilt to correctly 
represent the dielectric response of BZ at low tempera- 
ture. 

We now discuss the physical reasons behind the 
changes in the dielectric response induced by octahedral 
tilting. As noted above, the weaker dielectric response 
of the low-symmetry BaZrOa structure is caused by the 
40% drop in the contribution of the low frequency Last 



mode. To some extent, this effect can be rationalized 
by the rattling cation model. As the Og cages tilt, the 
O atoms are brought closer to the Ba atom. This re- 
duces the volume of the Ba-0i2 cage leaving less room 
for Ba to rattle around in response to an applied electric 
field. More precisely, the shorter Ba-O distances lead to 
a stiffening of the potential energy surface felt by the 
Ba. This gives rise to a smaller dielectric response, as 
is proportional to l/uP'. However, examination of Ba- 
O modes shows that lo shifts up by about 10%. This 
should lead to about 27% decrease in e^, accounting for 
about half of the Ba-O mode e^t change. The other factor 
that weakens the dielectric response is the decrease in the 
mode effective charge Z*, which is significantly smaller 
for the low-symmetry structure than for the 5-atom high- 
symmetry structure. Since Eqns. 3 and 4 shows that 
is proportional to the square of Z^^, distribution of Z* 
over multiple modes leads to a 20% decrease in e^. 

Comparison of the Ba-O phonon eigenvectors for the 
two structures shows that induced Zr off-centering is 
smaller for the Last mode in the 40-atom low symme- 
try structure than in the 5-atom high-symmetry mode. 
It is this reduction in the induced Zr displacement that 
gives rise to the smaller Z*. Physically, the smaller in- 
duced Zr displacement points to a weakened coupling be- 
tween Ba and Zr displacements. The coupling between 
Ba and Zr displacement is most likely due to Pauli re- 
pulsion |1^ . In the low-symmetry structure Og tilting as 
well as increased O motion in x and y directions in the 
Last mode eigenvector moves an O atom between Ba and 
Zr. This partially screens the Ba-Zr interactions, leading 
to a smaller Z* . 
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IV. CONCLUSIONS 

We have presented calculations that reveal interest- 
ing structural and dielectric properties in the seemingly 
simple material BaZrOa. At zero temperature, the high- 
symmetry model proved inadequate for calculating the 
dielectric constant. Large supercells were used in order 
to find the equilibrium structure, which exhibits octa- 
hedral tilts. These tilts have a significant effect on the 
low frequency Ba-0 Last mode and its contribution to 
the dielectric constant. Allowing octahdral tilts in the 
first principles modeling of BaZrOs resolves the discrep- 
ancy between recent theoretical and experimental inves- 
tigations of low temperature dielectric constant. More- 
over, the understanding of how octahedral tilitng can 
change dielectric response should give insight into alloy- 
ing strategies for future dielectric materials design. 
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